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ABSTRACT 

We report on combined INTEGRAL and XMM-Newton observations of the super- 
giant fast X-ray transient IGR J16418— 4532. The observations targeted the X-ray 
eclipse region of IGR J16418— 4532s orbit with continuous INTEGRAL observations 
across ^25% of orbital phase and two quasi-simultaneous XMM-Newton observations 
of length 20 ks and 14 ks, occurring during, and just after the eclipse respectively. 
An enhanced INTEGRAL emission history is provided with 19 previously unreported 
outbursts identified in the archival 18— 60keV data set. The XMM-Newton eclipse 
observation showed prominent Fe-emission and a flux of 2.8x 10 -13 ergcm" 2 s _1 (0.5 
- 10 kcV). Through the comparison of the detected eclipse and post eclipse flux, 
the supergiant mass loss rate through the stellar wind was determined as M w = 2.3- 
3.8x 10 -7 M yr _1 . The post eclipse XMM-Newton observation showed a dynamic flux 
evolution with signatures of the X-ray pulsation, a period of flaring activity, structured 
riff variations and the first ever detection of an X-ray intensity dip, or 'off-state', in a 
pulsating supergiant fast X-ray transient. Consideration is given to the origin of the 
X-ray dip and we conclude that the most applicable of the current theories of X-ray 
dip generation is that of a transition between Compton cooling dominated and radia- 
tive cooling dominated subsonic accretion regimes within the 'quasi-spherical' model 
of wind accretion. Under this interpretation, which requires additional confirmation, 
the neutron star in IGR J16418— 4532 possesses a magnetic field of ~10 14 G, providing 
tentative observational evidence of a highly magnetised neutron star in a supergiant 
fast X-ray transient for the first time. The implications of these results on the nature 
of IGR J16418-4532 itself and the wider SFXT class are discussed. 

Key words: X-rays: binaries - X-rays: individual: IGR J16418— 4532 - stars: winds, 
outflows - stars: pulsars: accretion, accretion discs 



1 INTRODUCTION 
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Supergiant Fast X-ray Transients (SFXT) are a sub-class 
of supergiant high mass X-ray binaries (HMXB) that dis- 
play extreme flaring behaviour on short (~ hour) timescales 
(Sguera et al. 20051. They also display an X-ray dynamic 
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range in excess of that possessed by classical supergiant X- 
ray binaries (~ 10-20, Sg-XRB), reaching 10 4 -10 5 in the 



most extreme systems such as IGR J17544— 2619 ( jRampy, 

Smith 
-i 



& Negueruela 2009) whilst only reaching about 
intermediate SFXTs, such as IGR 
Due to the supergiant 



10 z in the so-called 
J16465-4514 flClark et al.||2010l 



nature of their companion stars (Negueruela et al. 20061 



SFXTs are located along the Galactic Plane where there 



are currently 10 spectroscopically confirmed systems ( Sidoli 



2011 ), as well as a similar number of candidate SFXTs that 



show the required X-ray flaring behaviour but, as yet, do not 
have a spectroscopically confirmed supergiant counterpart. 

IGR J16418-4532 was first detected by INTEGRAL 
during observations of the transient black hole system 4U 
1630—47 (Tomsick et al. 20041 and proposed as a mem- 
ber of the SFXT class by ISguera et all ([20061 after the 



identification of short duration outbursts from the system. 



XMM-Newton observations reported by Walter et al. ( 2006 1 



identified a neutron star pulse period of 1246 ± 100 s (later 
refined to 1212 ±6 s by Sidoli et al. 20121 along with a 
high (~ 10 23 cm -2 ) absorption intrinsic to the system. The 
source was also localised to arcsecond accuracy and the IR 
counterpart proposed as 2MASS J16415078-4532253, which 



was later confirmed by Romano et al. ( 2012| through a 
Swift /XRT source localisation of RA,Dec (J2000) = I6 h 
41 m 50 s . 65, -45° 32' 27". 3 with an uncertainty of l".9 (90% 



confidence). Corbet et al. ( 2006[ ) identified a short, eclips- 
ing orbit of 3.7389±0.0004days and 3.753±0.004 days in the 
system using RXTE / ASM and Swift /BAT data sets respec- 
tive!^] (this determination was recently improved to 3.73886 
± 0.0002 8 days using RXTE / ASM data over an ~ 14 year 
baseline ( |Levine et al. 2011[ )). Combining the orbital and 
pulse periods places IGR J16418— 4532 in the wind-fed Sg- 
XRB r egion of the Corbet diagram ( |Corbet||l"986"| ). |Chaty| 
et al. ( 2008 1 performed NIR photometry on the IR coun- 



terpart of IGR J16418-4532 and, using SED fitting, de- 
termined a stellar temperature of 32,800 K suggesting the 
companion is of spectral type OB. Along with the location 
of IGR J16418-4532 in the Corbet diagram, the SED results 
confirmed the HMXB nature of this system and implied a 
minimum distance of ~ 13 kpc when assuming a supergiant 
companion. The supergiant nature of the companion was 
recently confirmed by Coleiro et al. (2013, private commu- 
nication) who classified the companion as a BN0.5Ia nitro- 
gen rich supergiant, thus confirming IGR J16418— 4532 as a 
intermediate SFXT. 

The maximal X-ray dynamic range observed in IGR 



J16418— 4532, greater than 10 in the soft X-ray regime ( Ro- 



mano et al. 2012 Sidoli et al. 20121, is in excess of that 



observed from classical, wind-fed Sg-XRBs with comparable 
orbital parameters. Through modelling the flare - luminosity 
distribution observed during an intensive Swift /XRT moni- 
toring campaign of IGR J16418-4532,|Romano et al.|f20l2l 



concluded that the enhanced level of variability in this sys- 
tem could be described by the NS accreting stellar wind 
clumps in the mass range 10 16 - 10 21 g embedded in a highly 
structured wind with a terminal velocity between 800 and 
1300 km s _1 (using the Ducci et al. ( 2009 1 formulation of the 



1 the statistical discrepancy in these two values is likely due to 
an underestimate of either one or both uncertainties 



'clumpy wind' scenario of SF XT outbursts ( |in't Zand|2005l 
Walter fc Zurita Heras|2007[ ) ). However using a 4 ks XMM- 
Newton observation in 2011, Sidoli et al. (20121 argue that 
the observed X-ray variability, orbital dynamics and quasi- 
periodic flaring activity is indicative of the NS accreting in 
a regime that is transitional between pure wind accretion 
and full Roche lobe overflow (RLO). It is argued that such 
a 'Transitional Roche Lobe Overflow (TRLO)' regime could 
be the dominant source of the observed X-ray variability in 
SFXTs with short orbital periods. 

In this work we present new quasi-simultaneous, or- 
bital phase targeted INTEGRAL and XMM-Newton obser- 
vations, along with an archival INTEGRAL study, of IGR 
J16418— 4532. In Section [2] we outline the data sets, analy- 
sis and results from both the archival and new INTEGRAL 
observations. Section [3] presents the analysis of, and results 
from, the new XMM-Newton data set. These results are then 
discussed in Section [4.41 and final conclusions drawn in Sec- 
tion [5] 



INTEGRAL DATA ANALYSIS AND 
RESULTS 

The archival data set from I NTEGRAL/IBIS ( | Winkler et aT] 



2003 Ubertini et al. 2003 1 consisted of all observations of 



IGR J16418-4532 spanning 2003 January 11 through 2010 
September 30 providing a total exposure of ~10Ms. All ob- 
servations were processed with version 9 of the INTEGRAL 



Offline Science Analysis software (OSA, Goldwurm et al 
12003 1 and images were created in the 18— 60keV energy 



range for each science window (ScW, individual INTEGRAL 
exposures of nominal length 2000 s). An 18— 60keV light 
curve was generated by extracting count rates and errors 
from each image at the best determined X-ray position of 
IGR J16418-4532. Our new INTEGRAL observations were 
performed between UTC 09:57:33 2012-09-01 and 09:00:29 
2012-09-02 for a total exposure of ~78ks. Images were again 
generated in the 18—60 keV band and a light curve extracted 
following the same procedures outlined for the archival data. 

The light curve was tested for periodicities by construct- 
ing a Lomb-Scargle periodogram ( |Lomb|1976| |Scargle|1982| 
and the peak period was identified as 3.74 days. The un- 
certainty on the peak was derived via a Monte-Carlo based 
randomisation test whereby the count rate of each point in 
the light curve is randomly varied within its error bar us- 
ing a Gaussian probability distribution. A periodogram is 
then generated for the altered light curve and the strongest 
period recorded. This process is repeated 50,000 times and 
the resulting distribution of peak periods fit with a Gaus- 
sian curve, the width of which is taken as the la error on 
the identified periodicity. The periodicity identified in the 
archival data set is therefore 3.7399 ± 0.0003 days. This 
value is consistent, within 2a, with the orbital period of 
3.73886 ± 0.00028days reported by |Levine eFaT] poTT| ), 
which was identified in an RXTE j ASM data set spanning 



14 years. As the Levine et al. (20111 determination has a 



slightly better constraint and originates from a longer data 
set it is this value of the orbital period that will be used for 
the remainder of this work. Additionally to maintain consis- 
tency with previous works the zero phase ephemeris of MJD 
53560.200 ( |Corbet et al.|2006[ > will be used throughout. 
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Figure 1. Archival INTEGRAL /IRIS 18 - 60keV light curve 
(scaled by a factor of 10) folded on P or b = 3.73886 days with the 
zero phase ephemeris MJD 53560.20. The new IBIS observations 
are shown by the green points and the phase locations of the 
XMM-Newton observations are shown in red. 



Figure 2. Orbital Phase Distribution of the archival IGR 
J16418— 4532 outbursts using the ephemeris P or b = 3.73886 days 
and a zero phase of MJD 53560.20. 



Figure[T]shows the archival 18 — 60 keV IBIS light curve 
folded on the 3.73886 day orbital period and scaled by a 
factor of 10 for illustrative purposes. Additionally the 18 - 
60 keV ScW light curve of the new INTEGRAL observations 
is over plotted in green and the orbital phase locations of 
the XMM-Newton observations are shown in red (see below 
for further details). The orbital profile is dominated by the 
deep eclipse that spans ~0.2 of the phase space between <j> 
= 0.45 and 0.65, corresponding to a duration of ~0. 75 days, 
and is consistent with the neutron star being fully eclipsed. 
Outside of the eclipse, the orbital profile is consistent with a 
constant flux at an unsealed IBIS count rate of ~1. 2 counts 
per second, corresponding to a flux of 7 mCrab in the 18 — 
60 keV band, showing that the eclipse is the driving factor 
in the observed periodicity. 

The archival 18 — 60keV light curve was also searched 
for outbursts by systematically identifying regions of the 
light curve with a high local significance in windows of 0.02 
to 3.0 days (~half an hour to the majority of one orbit). 
Due to the coded aperture nature of IBIS and the deconvo- 
lution method of image reconstruction, light curves of tran- 
sient sources contain a high amount of systematic noise that 
is symmetric about zero during times of non-detection. To 
guard against false outburst detections resulting from the 
large number of trials incurred by this method of outburst 
identification, the count rates in the light curve are inverted 
about zero and the same identification procedure performed. 
The maximum significance detected in the inverted light 
curve, which is a false positive by definition, is then taken 
as the cut-off point below which the local significances can 
not be safely distinguished from random statistical varia- 
tions present in the light curve. Applying this cut-off to the 
original set of excesses produced 35 distinct outburst events 
within the archival light curve down to a minimal signifi- 
cance of 4.48ct and with durations in the range ~0.02 to 
1.3 days. 16 of the 35 detected outbursts have already been 
reported in previous works (see |Sguera et~a l. 2006 , Ducci, 
Sidoli, & Paizis 2010 1 while the remaining 19 are new de- 



tections, the main properties of which are reported in Table 

m 

The total duration of these outbursts places the ac- 
tive duty cycle of IGR J16418-4532 at 6.14%. This figure 
is somewhat higher than that derived from the results of 



Ducci, Sidoli, & Paizis (20101 which give a duty cycle of 



-4.3%. Whilst the analyses were performed on data from 
different energy bands (18 — 60 keV here and 20 — 40 keV in 



Ducci, Sidoli, & Paizis 2010 1 and on data sets of different 



total exposure, it is likely the difference in the methods used 
to identify outbursts in each work is the main cause of this 



disparity. Ducci, Sidoli, & Paizis (20101 consider only the 



ScWs in which IGR J16418— 4532 was detected at a signif- 
icance of > 5cr from which to define outbursts whereas here 
we consider all ScWs in our searches and model the outburst 
significance cutoff to apply as described above. Under this 
method an outburst search becomes sensitive to a popula- 
tion of longer, fainter outbursts which are below the formal 
detection threshold for a known source in a single ScW (i.e. 
5a) but sum to a significant detection when the multiple 
ScWs of an outburst are combined. The extent to which 
the size of the increase in the active duty cycle observed for 
IGR J16418— 4532 is specific to this source or represents a 
more global phenomenon however requires further investi- 
gation. Figure [2] shows the distribution, in orbital phase, of 
the peak emission from each identified outburst. It can be 
seen that, apart from the eclipse region where no outbursts 
are detected, there is no coherent structure in the outburst 
distribution which is consistent with the phase folded light 
curve (Fig. [I]). This implies that outbursts are equally likely 
to occur at any point in the orbital phase. 

The new INTEGRAL observations were constrained 
in orbital phase to cover the eclipse region of IGR 
J16418-4532's orbit, covering <j> = 0.403 through 0.655 us- 
ing the above ephemeris. The light curve is shown by the 
green points in Fig. [T] and it can be seen to cover the eclipse 
ingress and fully eclipsed region but ends just before the 
egress. IGR J16418— 4532 shows a low, but rising level of 
activity in the first 4 ScWs before the eclipse which is then 
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56171.55 
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Figure 3. 1000 s binned 18 — 60 keV light curve of the low level 
activity observed during the first 5 ScWs of the recent INTE- 
GRAL/IBIS observations. 



cutoff sharpfy in the fifth ScW at an orbital phase consistent 
with the eclipse ingress. The source is not detected during 
the subsequent ScWs performed during the eclipse. A mo- 
saic image of the first 5 ScWs was produced in which IGR 
J16418— 4532 was detected at a significance of 6.1a. A spec- 
trum was extracted from these ScWs and fit with a powerlaw 



in XSPEC version 12.7.1 (|Arnaud||l996|) with V = 2.2+ 1 1 7 



and x =1.08(6 do}). The correspoITHTng 18-60keV flux 
was 1.4 x 10 -10 erg cm -2 s _1 . A higher time resolution light 
curve with 1000 s bins was also extracted from these ScWs. 
This light curve is shown in Fig. [3] and it again shows a 
low level of activity, followed by an increase in flux at MJD 
56171.53 which then rapidly drops again. The peak emis- 
sion occurs at MJD 56171. 537±0. 0058, relating to an orbital 
phase of <j) = 0.431. Whilst it is possible that the flux varia- 
tion could be due to X-ray flaring the decrease in flux could 
also be resulting from eclipse ingress which would place a 
limit on the ingress time of ~ 2000 s (Acf> = 0.006) and pro- 
vides an accurate eclipse ingress ephemeride. Unfortunately 
the signal-to-noise was insufficient to probe this emission at 
higher time resolutions. 

IGR J16418-4532 was not detected in the soft 
3 -lOkeV X-ray band by the JEM-X instrument ( |Lund et| 
al. 2003 1 aboard INTEGRAL during these recent observa- 



tions. The 3a upper limit derived from the mosaic images 
of both JEM-X units for an exposure time of 14.5 ks, the 
majority of which occurred whilst IGR J16418— 4532 was in 
eclipse, is 3.4mCrab (3— lOkeV). 



XMM-NEWTON DATA ANALYSIS AND 
RESULTS 



Two separate XMM-Newton /EPIC (Jansen et al. 2001 



Turner et al. 2001 Striider et al 



2001) observations 



were performed quasi-simultaneously with INTEGRAL to 
achieve sensitive soft X-ray coverage over a wider region 
of orbital phase. The first observation was performed be- 
tween UTC 21:12:25 2012-09-01 and 02:53:22 2012-09-02 for 
an exposure of ~20ks. This observation covered the orbital 



phase range <j> = 0.525 — 0.587 and occurred whilst the NS 
was fully eclipsed by the supergiant. The second observation 
was performed after the NS had egressed from the eclipse, 
between UTC 14:16:33 and 19:12:03 2012-09-02, for an ex- 
posure of ~14ks and covered the orbital phase range cf> — 
0.715 — 0.758. The phase location of these observations is 
shown in Fig. [I] In both observations the EPIC-MOS and 
EPIC-pn detectors were operating in large window mode. 

Data from both the EPIC-MOS and EPIC-pn detec- 
tors were analysed from each observation using SAS vl2.0.1 



(Gabriel et al. 2004) and the most recent instrument cali- 



bration files. The data sets were checked for regions of high 
particle background following the method outlined in the 
XMM-Newton SAS data analysis thread^ with no regions 
of high particle background identified in either the MOS or 
pn 10 — 12 keV light curves using cut-offs of 0.35 and 0.4 
counts s' 1 respectively. The EREGIONANALYSE tool was 
used to define the optimal extraction region for all light 
curve and spectral generation procedures (note: all subse- 
quent references to optimal extraction regions were defined 
using this method). The data sets were evaluated for the 
presence of photon pile-up using the EPATPLOT tool and 
by comparing the shape of spectra extracted from both cir- 
cular and annular extraction regions. The observations were 
found to be unaffected by pile-up apart from the final few 
ks of the second data set, during which the shape of the 
extracted spectra was observed to vary slightly when circu- 
lar and annular extraction regions were used. As such the 
standard annular extraction region method was utilised to 
remove the effects of pile-up from the extracted X-ray spec- 
tra during these times. The specific extraction region shapes 
and sizes used for both data sets are outlined in Sections l3.1 
and |3.2| All spectra reported in this work were extracted 
following the standard procedures and the SAS tools RM- 
FGEN and ARFGEN were used to extracted the necessary 
response files for each spectrum. Spectra were again fit us- 
ing XSPEC version 12.7.1 with uncertainties quoted at the 
90% confidence level throughout and the elemental abun- 



dances set to those of Wilms, Allen, & McCray (20001. Due 



to the very different nature of the detected emission in the 
two observations, the detailed analysis of each is reported 
separately in Sections |3. 1| and |3.2| below. 



3.1 Mid-eclipse Observation 

Emission from IGR J16418— 4532 was detected at a low level 
in the mid-eclipse observation (MEO) which, given its loca- 
tion in orbital phase, is assumed to originate from repro- 
cessing of the the NS X-ray flux by the stellar wind of the 
supergiant. There was insufficient signal to produce mean- 
ingful light curves and spectra in the MOS detectors so they 
will not be considered further. Similarly in the pn detector 
only a low signal-to-noise light curve could be extracted in 
a broad (0.2 — lOkeV) energy band. The light curve did 
not provide any detailed temporal information, however a 
spectrum was accumulated from the full exposure of the pn 
detector using an optimal circular extraction region of ra- 
dius 14". The background subtracted pn spectrum of the 
MEO, binned to a minimum of 25 counts per bin, is shown 



http: / / xmm.esac.esa.int /sas / current /documentation/threads / 
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Table 1. Newly discovered outbursts of IGR J16418— 4532 in the archival INTEGRAL/IBIS data set. The orbital phase given relates to 
the ScW with the maximal count rate in each event. The peak flux is that of the same ScW in the 18-60keV band, calculated assuming 
a powerlaw spectrum with T = 2.2 as observed in the activity at the start of the new INTEGRAL/IBIS data set. Note: Fluxes denoted 
by a * are likely over estimates of the source flux due to large statistical uncertainties in the detected count rate. 
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in the top panel of Fig. [4] It can be seen that the spectrum 
has a flat shape with a large emission feature consistent with 
a 6.4 keV iron-Ka line. 

To further characterise the spectral shape of the pn 
spectrum it was first fit with a simple absorbed power- 
law continuum with an additional Gaussian component 
(PHABS(POWERLAW + GAUSSIAN) in XSPEC). The 
absorption was fixed to the Galactic v alue in the direction 
of IGR J16418-4532, 1.88xl0 22 cm -2 (|Dickey fc Lockman 



19901, as in this observation the powerlaw is simply param- 
eterising a more complicated reprocessing mechanism and is 
unlikely to be describing the intrinsic, locally absorbed con- 
tinuum emission. The power law index, Gaussian line energy 
and line width were left as free parameters. The spectrum 
was well fit by this model with \ 2 = 1.02 (8 dof), however 
the Gaussian line width o~e = 0.15 keV (EW =3.1 keV) cor- 
responds to an electron plasma temperature of ~3xl0 6 K 
which is well in excess of that expected for a supergiant 
atmosphere. This unphysical fit is likely resulting from the 
low signal-to-noise and coarsely sampled spectrum combined 
with the possible influence of the 7.1 keV iron absorption 
edge located within the emission line region of the spectrum. 

To obtain a more physical characterisation of the iron 
line profile the width of the Gaussian line was fixed at 
such that the line is intrinsically narrow and broadened only 
by the EPIC-pn response. Using a single Gaussian with the 
centroid energy as a free parameter provided a poor fit to the 
data (x 2 = 2.04 (9 dof)) with large residuals in the region 
of the emission line. A second intrinsically narrow Gaussian 
was added at higher energy to account for emission from 
ionised iron and the spectrum fit with the centroid energy 
left as a free parameter. The resulting fit was of good quality 
with \ 2 = 0.72 (7 do}) and the line centroid energies Ei 
= 6.40 t°QM and E 2 = 6.65 t°o° t keV, consistent with both 
neutral 6.4 keV and ionised 6.67 keV Fe-Ka emission. The 



Table 2. IGR J16418-4532 MEO best fit EPIC-pn spectral pa- 
rameters and uncertainties. The spectrum was fit using the model 
PHABS(POWERLAW + GAUSS + GAUSS) with n H fixed and 
X 2 = 0.72 (7 dof). The flux is quoted in the 0.5-10keV band 
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equivalent widths of the 6.40 and 6.65 keV lines were 0.63 
and 0.33 keV respectively showing that the neutral emission 
is approximately twice the strength of the ionised emission. 
The 0.5 - lOkeV detected flux was (2.80 tolt) x 10~ 13 erg 
cm -2 s _1 . The best fit parameters resulting from this model 
are outlined in full in Table [2] with the best fit model and 
residuals shown in Fig. [4] where the lower panels illustrate 
the improvement in the residuals as the Gaussian lines are 
added to the model. 



3.2 Post-eclipse Observation 



During the post-eclipse 
J16418-4532 was detected 
Figure [5] shows the 0.2 — 



observation (PEO) IGR 
with good signal-to-noise. 
lOkeV EPIC-MOS1 (black) 
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Figure 4. Top panel: EPIC-pn best fit spectrum of IGR 
J16418-4532 during the MEO. The spectrum is best fit with an 
absorbed powerlaw with two intrinsically narrow Gaussian emis- 
sion lines at energies of 6.40 ^q'q4 an d 6.65 1.q'q5 keV respectively. 
Lower panels: The three lower panels show the residuals of the fit 
when using an absorbed powerlaw, an absorbed powerlaw with 
one intrinsically narrow Gaussian and the best fit absorbed pow- 
erlaw with two intrinsically narrow Gaussians. 



and MOS2 (red) light curves, extracted from optimally 
defined circular regions of radius 97" and 98" respectively, 
with time binning of 100 s. T zero is denned as the first 
time stamp in the MOS1 light curve and takes the value 
Tzero = MJD 56172.596. The light curves reveal a range 
of activity including a steady emission at ~ 2 counts 
s _1 , a rapid 'dipping' feature between t ~ 800 and 2500 s 
and flaring behaviour at times of t > 10,000 s. There is 
also activity on short, hundreds of seconds timescales 
throughout the exposure, some of which may be due to 
the ~1200s pulse period. This activity was also observed 
in the EPIC-pn detector, the 100 s binned light curve of 
which is shown in the top panel of Fig. [6] The pn light 
curve was extracted from an optimal circular region of 
radius 68". Due to additional observational overheads 
the pn detector did not start exposing until t = 1700 s, 
resulting in a lack of coverage during the majority of the 
'dip' region observed by the MOS detectors. To account 
for the effects of pile-up at times greater than t = 10,000 s, 
annular extraction regions with inner radii of 7.5" and 10" 



for the MOS and pn detectors respectively were used to 
extract all data products. 

The maximal dynamic range inferred from the 100 s 
binned broad band 0.2-10 keV light curves is 75 (MOS1) 
indicating a large variation in the flux detected during the 
relatively short exposure. 100 s binned light curves were also 
extracted in restricted 0.2—4 and 4— lOkeV energy bands 
and the 0.2—4 / 4— lOkeV hardness ratio (HR) calculated 
for each detector. The EPIC-pn HR is shown in the lower 
panel of Fig. [6] and it is seen to be composed of features 
on varying time scales. The most striking feature is the sys- 
tematic hardening of the detected emission between t = 6000 
and 10,000 s and the strong correlation between the HR and 
detected source flux minima. The end of this hardening is 
also temporally coincident with the onset of the flaring be- 
haviour observed from t = 10,000 s onwards. Variation is also 
observed on timescales of a few 100 s throughout the ob- 
servation, regardless of flux level, and there is evidence of 
systematic softening of the detected emission as the source 
exits the 'dip' region of the light curve. HR light curves were 
also generated for the MOS detectors that showed the same 
evolution through the observation, albeit with less clarity 
due to the larger count rate errors, but could not provide 
additional diagnostics during the 'dip' region as a result of 
the limited photon statistics at this time. 



3.2.1 Spectral Analysts 

The lower panel of Fig. [6] illustrates the level and speed of 
variation in the emission of IGR J16418— 4532 during the 
PEO through the evolution of the hardness ratio, imply- 
ing fast variations in the detected source spectrum. In order 
to track these variations we extracted spectra accumulated 
from bins of length 400 s in steps of 200 s through the length 
of the exposure. An image was extracted for each time bin 
and the optimal extraction region for that region defined 
by the combination of EREGIONANALYSE and the cir- 
cle/annulus shape constraints described above. Spectra were 
extracted, and binned to a minimum of 15 counts per bin, 
from both the MOS and pn detectors at times greater than 
t = 1700 s and the MOS detectors only prior to this. Such a 
systematic approach was adopted to remove possible biases 
introduced by the subjective selection of 'similar' regions 
while the extraction bin size of 400 s was chosen as a com- 
promise between generating spectra of sufficient signal-to- 
noise to perform accurate fitting and reducing the amount 
of HR variation within each bin. The spectra were binned at 
15 counts per bin to ensure a sufficient number of spectral 
bins to allow model fitting were present even in the spectra 
accumulated from the time periods of lowest flux. Step sizes 
of 200 s were chosen so that the evolution of spectral param- 
eters could be followed at sufficient resolution but the data 
would only be over sampled by a factor of 2. 

Inspection of the spectra showed no evidence of Fe 
emission lines so the spectra were fit with simple absorbed 
models to characterise their shapes. Absorbed powerlaws 
(PHABS(POWERLAW) in XSPEC) were simultaneously fit 
to the spectra from all available detectors in each time bin 
in the 0.5 — 15 keV energy range with the absorption, pow- 
erlaw index and normalisation left as free parameters. The 
evolution of the spectral parameters, uncertainties and de- 
rived fluxes across the PEO are shown in Fig. [7] The inset 
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Figure 5. EPIC-MOS1 (black) and MOS2 (red) 100 s binned light curve of IGR J16418-4514 in the 0.2 - lOkeV band 



panel in Fig. [7] shows the x 2 value of the fit to each bin 
throughout the observation. For the majority of the expo- 
sure time the x 2 values are distributed about 1.0 showing 
a good fit to the spectra. However bins falling within the 
'dip' region of the observation, in particular those occurring 
during the deepest part of the 'dip', gave unacceptable x 2 
values, indicating that an absorbed powerlaw is a poor de- 
scription of the emission at these times. The spectra of this 
region are considered separately and in more detail below. 
For the rest of the observation it can be seen that when both 
xih and F are left as free parameters there is a level of varia- 
tion in both parameters. However a Chi-Squared test to the 
mean value of each parameter, where both the mean and 
Chi-squared statistic have been calculated with the poorly 
fitted bins omitted, gave a x 2 value of 0.65 (67 dof) and 3.52 
(67 dof) for F and xih respectively. Hence the F parameter is 
consistent with being constant throughout the observation 
whilst the absorption shows excess variation. The average 
r value is 1.195 and is over plotted as the red line in the 
second panel of Fig. [7] 

A second set of fits was then performed with F fixed at 
the value of 1.195, whilst the normalisation and xxh were left 
as free parameters. The evolution of the spectral parameters 
and the goodness-of-fit in this case are shown in Fig. [8] A 
similar, better constrained xih evolution as that seen in Fig. 
[7] is again observed when fitting a powerlaw of fixed V with 
this model providing an acceptable fit for all bins apart from 
those occurring during the 'dip' feature. The xih evolution 
shows that there is a constant, high level of absorption of 
n H ~7xl0 22 cm 2 that is intrinsic to the IGR J16418-4514 
system (the Galactic value in the direction of the source is 
1.88x 10 22 cm -2 ). The most striking feature of the evolution 
is the large, structured increase in xih between t = 6500 and 



9500 s which is co-incident with the 'Pre-flare hardening' fea- 
ture observed in the HR. In particular the peak n^ is again 
temporally coincident with the HR and detected flux min- 
ima. Figure [9] shows an expanded view of the xih evolution 
overlaid with the unabsorbed flux from each bin, calculated 
assuming the fixed powerlaw of F = 1.195. It is seen that re- 
moving the effects of absorption brings the flux back up to 
the level seen before the nn increase but it is not until after 
the xih has returned to its base level that the large increase 
in flux starts to occur. Consideration of the physical inter- 
pretation of this behaviour is given in Section [4.4| Outside 
of the large increase the xih also shows some level of varia- 
tion that occurs on the same time scale as the NS pulsation 
(~ 1200 s). This variation is more pronounced in the spectral 
fits performed with a fixed F, adding to the justification for 
fixing the photon index as this produces spectral variations 
that may be linked to a physical origin whereas fixing the 
xxh and allowing Y to vary only produces fits with very poor 
X 2 values. 

The time resolved spectra were also fit with the more 
complicated 'Partial Covering' absorption model (PCFABS 
in XSPEC) in an attempt to better understand the evolution 
of the absorbing material during the PEO. However it was 
seen that the fits were insensitive to the additional parame- 
ters, with the covering fraction tending towards 100% for all 
spectra while the total xih (the combination of a component 
fixed to the galactic value and a free varying, partial cover- 
ing component) was consistent with the best fit value for the 
simple absorbed powerlaw model. Hence no additional infor- 
mation could be drawn from this model. The time resolved 
spectral analysis was performed with different bin durations, 
step sizes and over sampling factors and the evolution was 
observed to follow a consistent path regardless of the extrac- 
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Figure 6. Top: EPIC-pn 100 s binned light curve of IGR 
J16418-4532 in the 0.2-10 keV band. Bottom: EPIC-pn 0.2-4 / 
4— lOkeV hardness ratio light curve at the same resolution where 
lower values indicate a hardening of the detected flux 



tion used. As a final test a blackbody continuum was used in 
place of the powerlaw model. A consistent riu and flux evo- 
lution, with a similar distribution of \ 2 values, was observed 
throughout the observation producing further evidence that 
the evolution is explained by a relatively constant continuum 
modified by a varying photoelectric absorption. For the re- 
mainder of this work we refer only to the simple absorbed 
powerlaw fits when discussing the spectral evolution. 

To characterise the emission during the dip region of 
the light curve a separate EPIC-pn spectrum was accu- 
mulated encompassing all available EPIC-pn data of the 
region, namely between t = 1678 s and t = 2000 s for a to- 
tal exposure of 322 s. The spectrum showed no indica- 
tion of Fe emission lines and was fit with a simple ab- 
sorbed powerlaw where V was fixed to 1.195 (as above) and 



the best fit n H = (11.1] 



xW 2 



The x 2 = 0.72(15 



dof) and the corresponding 0.5— lOkeV unabsorbed flux 
was 2.0xl0 -11 ergcm -2 s _1 . When compared to the max- 
imum, unabsorbed flux detected later in the observation 
(27.2xl0" 11 ergcm- 2 s- 1 (0.5 - lOkeV)) a soft X-ray dy- 
namic range of 14 is observed during the PEO. 



3.2.2 Timing Analysis 

To test for signatures of the known 1212±6s NS spin period 
( Sidoli et al.|2012[ ) the 0.2-10 keV light curves from each in- 
strument were re-extracted with a variety of bin sizes. Here 
we report on the results from light curves with a binning 
of 10 s but point out that the results were consistent in the 
analysis of each different binning. For the purposes of tempo- 
ral analysis the 'dip' region of each light curve was removed 
as it occurs over a similar timescale to the known pulsation 
and could therefore generate significant power at frequen- 
cies similar to that of the NS pulse period, despite it not 
necessarily being of a rotational nature. The Lomb-Scargle 
technique was again applied to the light curves and strong 
signals detected at ~403 s in all detectors. From here we con- 
sider only the MOS2 light curve as it has a longer baseline 
than that of the pn and the results are consistent with those 
from MOS1. The uncertainty on the peak was calculated 
using the method outlined in Section [2] and the peak period 
found to be 403.04±0.14s. This is consistent with being the 
third harmonic of the known pulse period with the funda- 
mental frequency relating to a period of 1209.12±0.42 s. This 
spin period is consistent with the previously reported values, 



summarised in Sidoli et al. (2012), and with the addition of 



this detection there is still no significant detection of spin 
period evolution in IGR J16418-4532. 

The upper and middle panels of Fig.[l0]shows the MOS2 
light curve folded on the fundamental period of 1209.12 s in 
the 0.2— 4keV and 4— lOkeV bands respectively. A similar 
pulse profile, with one sharp narrow peak and one equally 
high but broader peak, is detected in both energy bands 
with the 0.2— 4keV displaying a pulse fraction of 25.0±3.6% 
and the 4-10 keV band 37.2±2.7%. The pulse fraction of the 
folded profile in the broad 0.2-10 keV band was 31.2±2.1%. 
A double peaked profile is consistent with the 2011 XMM- 



Newton observations of Sidoli et al. (2012 I and shows a con- 



tinued departure from the 2004 XMM-Newton observations 



of | Walter et al. ( 2006|) which sho wed a pulse profile with a 
single, broad peak ( |Sidoli et aT 2012 1 . The double peaked 



profile with uneven pulse widths could be the driving factor 
behind the detection of the third harmonic, as opposed to 
the fundamental frequency in these observations. The lower 
panel of Fig. [To] shows the hardness ratio as a function of 
pulse phase and it is seen that there is no systematic varia- 
tion of the hardness with the progression of the pulsation. 



4 DISCUSSION 

4.1 Stellar and Orbital Parameters 

IGR J16418-4532 is an intermediate SFXT that has shown 
some peculiarities in its emission history. Here we have pre- 
sented new combined INTEGRAL and XMM-Newton obser- 
vations of the source that targeted the eclipse region of IGR 
J16418— 4532s orbit. The well determined orbital period of 
3.73886±0.00028days ( |Levine et al.|2011[ ) allows us to place 
constraints on the orbital geometry and nature of the com- 
panion star. The flat, unstructured phase-folded light curve 
and orbital phase distribution of outbursts (outside of the 
X-ray eclipse) detected by INTEGRAL, as shown in Fig. [I] 
and [2] respectively, indicate that the orbit is likely circular 
(e ~ 0) due to the lack of coherent modulation of the hard 
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Figure 7. Temporal evolution of the best fit spectral parame- Figure 8. Temporal evolution of the best fit spectral parameters 

ters of an absorbed powerlaw to the 400 s binned spectra. Top of an absorbed powerlaw with a fixed T = 1.195. The panels show 

to bottom: EPIC-MOS2 0.2-10keV light curve; Powerlaw in- the same plots as Fig. [7] 

dex evolution; njj evolution; inferred 0.5 — lOkeV absorbed flux 
where the diamond points represent fluxes derived from the spec- 
tra accumulated from annular extraction regions. Inset: Reduced 
X 2 values for the fit to each spectrum. 
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X-ray flux across the orbital phase. A similar flat orbital 
profile is also seen in IGR J16479— 4514 which has a short 
3.32 day eclipsing orbit (Jain, Paul, & Dutta |2009|). Con- 



versely coherent orbital phase modulation has been observed 
in other SFXT systems which show evidence of a significant 
eccentricity in their orbit, such as a strongly peaked or- 
bital emission profile with outbursts clustered around times 
of presumed periastron, i.e. IGR .117354—3255 (Sguera et 
aL||2011 1. Additional dynamical considerations also suggest 
IGR J 16418— 4532 should possess a circular geometry as 
with such a short orbital period the system would be ex- 
pected to circularise relatively early in the lifetime of the 
X-ray binary. 

The stellar parameters of the supergiant companion in 
IGR J16418— 4532 have very few dynamic or spectroscopic 
constraints however. By combining the spectral classifica- 
tion (Coleiro et al. 2013, private communication), the IR 
SED fit parameters ( Chaty et al.|2008 l, the well determined 
orbital period ( jLevine et al. 2011[) a nd the the LI point sepa- 
ration equations of |Eggleton| (|1 983| we can attempt to place 
some constraints on the orbital parameters through the con- 
dition that full RLO does not occur in IGR J16418-4532. 
Throughout these considerations we assume a neutron star 
in a circular orbit with a mass of 1.4 Mq, a radius of 10 km 
and a source distance of 13kpc (Chaty et al. 20081. We 



utilise the Rs g /Ds g SED parameter of |Chaty et al.| ( |2008[ ) to 
calculate the range of possible stellar radii of the companion 
star, the central value of which (3.77 x 10~ 1 ) gives a stellar 
radius of 21.7 R0. However, calculations of the LI point sep- 
aration for a range of companion star masses indicate that 
RLO would be occurring for all companion masses of less 
than 50 Mq. Alternatively the lower limit of the Rsg/Dsg 
parameter (2.64 xlO -11 , 90% confidence) gives a stellar ra- 
dius of 15.2 R© and RLO is only initiated up to a stellar mass 
of 20 M©. For the remainder of this discussion we assume a 
stellar mass of 30 M©, the LI point separation of which is 
17.8 R©, and a stellar radius of 17 Rq such that the super- 
giant almost fills its Roche Lobe. Under these parameters 
the orbit of IGR J16418— 4532 is characterised as circular, 




Figure 10. Pulse phase folded profiles using the pulse period 
of 1209.12 s where the zero phase is arbitrary. Top: 0.2— 4keV 
folded profile. Middle: 4-10keV folded profile. Bottom: HR of 
folded profile. 



at an approximately edge on inclination with a radius of 
31.9R0 and a NS orbital velocity of 432 kms -1 . 

While the choice of these parameters may seen some- 
what overly arbitrary we wish to point out that the defining 
parameter in our remaining considerations, due to its strong 
influence on the size of the accretion radius, is the velocity 
of the stellar wind at the orbital separation of the neutron 
star. Following previous works on the properties of the stel- 



lar wind in IGR J16418-453 2 (i.e. |Romano et all pQ12| |, 
Ducci, Sidoli, fc Paizis| (|2010[)) we adopt the velocity law of 
Castor, Abbott, & Kleinl (119751) to describe the stellar wind 



velocity for both the smooth and clumped wind components, 
namely 



v(r) 



0.9983 



Rs 



9^ 



(1) 



The stellar mass influences this relationship through the de- 
rived value of the orbital separation r, and the stellar radius 
influences through the value of Rob ■ However equally influ- 
ential are the values of the terminal wind velocity Voo and 
the powerlaw index /?. Through modelling of the flare lu- 
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minos i ty di stribution as observed by Swift /XRT, Romano 



et al. 



(20121 place Voo in the range 800 to 1300 km s -1 and 
(3 between 0.8 and 1.3. The value of the stellar wind veloc- 
ity at the orbital separation calculated with the above stel- 
lar parameters and «oo = 1000 km s" 1 , P = l is 470 km s -1 . 
However by varying the «oo and /3 parameters within their 
limits the stellar wind velocity occupies the range ~370 to 
610kms _1 . Similarly if we select a different set of allowed 
stellar parameters (Ms 9 = 4OM0, Rs 9 = 2OR0) and calcu- 
late the stellar wind velocity for the same range of allowed 
Voo and ft values the wind velocity at the orbital separation 
is in the range ~340 to 560 km s -1 . Hence the largest source 
of uncertainty in the stellar wind velocity comes from the 
poorly constrained wind parameters and not the unknown 
stellar parameters. For the remainder of this discussion we 
therefore adopt the stellar parameters outlined above and 
calculate the relative velocity of the NS orbital and stellar 
wind radial velocities as v re i = (v% s +«^)5. Uncertainties 
in derived physical values are then discussed in terms of the 
range of possible stellar wind parameters. 



4.2 Eclipse Emission 

Assuming that the cut off in flux observed by INTE- 
GRAL/IBIS was due to the eclipse ingress this provides an 
accurate ephemeris of MJD 56171. 537±0. 0058 for the start 
of the X-ray eclipse. Combining this with the first XMM- 
Newton observations (MEO) allows a firm lower limit, of 
0.583 days (A0 = 0.156), to be placed on the duration of 
the X-ray eclipse in IGR J16418-4532. 

The detection of X-ray flux from HMXB systems whilst 
the X-ray source is eclipsed by the supergiant companion 
is a common phenomenon amongst such sources, being de- 



tected in both classical Sg-XRBs (e.g. Cen X-3, |Naik, Paul, 
fc Ali|201lD and SFXTs (e.g. IGR J16479-4514, [Bozzq et aL 
2008[). The detected 0.5-10 keV flux of 2.8x10 "erg cm -2 



s~ observed during the MEO corresponds to a luminosity 
of 6.3xl0 33 erg s _1 at the assumed source distance, ruling 
out direct emission from the companion star as the source of 
the detected emission. Instead the detections during eclipse 
are usually interpreted as resulting from the reprocessing of 
the intrinsic flux from the NS through the dense wind of 
the supergiant. In the MEO emission lines at energies con- 
sistent with atomic transitions from both neutral (6.4 keV) 
and ionised (6.67 keV) Fe were detected. The lack of Fe line 
features in the PEO indicate that the iron emission does 
not originate from close to the NS, but in extended re- 
gions throughout the supergiant atmosphere, with the rela- 
tive equivalent widths of the lines (0.63 and 0.33 keV for the 
6.4 and 6.67 keV lines respectively) suggesting that neutral 
Fe is the more abundant state. 



4.3 Out of Eclipse Emission 



Following the method of Sidoli et al. (20131, it is possible 



to use both XMM-Newton observations (MEO and PEO) 
to probe the supergiant total mass loss rate and the stel- 
lar wind density at the orbital separation. Using the system 
parameters discussed above gives an orbital separation of 
2.2xl0 12 cm. Additionally the average flux detected whilst 
the NS was eclipsed in the MEO (2.8x 10~ 13 erg cm" 2 s _1 ) is 



~0.56% of the average PEO flux detected outside of the X- 
ray dip and flare regions (~ 5. Ox 10 -11 erg cm -2 s" 1 ). This 
allows the wind density to be estimated, assuming that the 
X-ray emission detected during the eclipse is produced by 
X-ray scattering, as n w (a) = 0.0056/(aoT), where <tt is the 
Thomson scattering cross section and a is the orbital sepa- 



ration (Lewis et al. 1992) 



The wind density at the orbital separation, p w (a), is cal- 
culated as 6.3x 10 -15 gem -3 . Assuming that the wind has 
a spherical geometry, and 13=1 (Castor, Abbott, & Klein 



19751, the ratio between the wind mass-loss rate and the 



terminal velocity can be calculated as M w /v oc =4 : ty a (a — 
R op t)p^(a), giving A^/u oo =9.3xlO- 18 M km- 1 . If the 
wind terminal velocity is in the range 800-1300 km s" 1 (fRa 



mano et al.|2012| |, the mass loss rate via the stellar wind lies 
in the range M w = (2.3-3.8) x 10~ 7 M yr" 1 . Assuming di- 
rect accretion from the wind, the accretion rate onto the NS 
can be estimated as M acc = ((i"^ 2 cc /47ra 2 )M 11 ,), where R acc 
is the accretion radius and is given by 2GMjvs/v 2 ei . This 
implies an accretion rate at the orbital separation of be- 
tween M acc ~ (9.7 - 5.4) xlO 15 gs _1 , corresponding to X-ray 
luminosities in the range, Lx = (1.8— 1.0) x 10 36 ergs -1 . 

The average, unabsorbed flux detected during the PEO, 
outside of the X-ray dip and flare, corresponds to a lumi- 
nosity of 2. Ox 10 36 ergs -1 at a distance of 13kpc, imply- 
ing that the estimated wind density at the orbital separa- 
tion, for the range of parameters assumed here, is able to 
generate the accretion luminosity observed outside of the 
eclipse in IGR J16418— 4532. This is in contrast to the case 
of IGR J16479— 4514 where the calculated stellar wind den- 
sities produced an X-ray luminosity two orders of magnitude 
greater than that observed, suggesting the need for a damp- 
ing mechanism to reduce the effective mass accretion rate 
in this system ( Sidoli et al.|2013 l. Such a mechanism would 
appear to not be required at nominal luminosities in the case 
of IGR J16418-4532. 



4-3.1 tih variations 

The evolution of the spectral parameters and source flux ob- 
served during the PEO (Fig. [sj is intriguing, with evidence 
of a very variable accretion environment being experienced 
by the NS in IGR J16418-4532. Two striking features are 
observed during the course of the observation, the first be- 
ing the deep dip observed in the detected flux that occurs 
between ~ t = 800 and 2500 s which is discussed separately 
below. The second is the pronounced increase and evolution 
in absorption observed between t~ 7000 s and t~ 10,000 s. 
The absorption first shows a sharp rise to ~ 10 23 cm -2 from 
the base level of ~ 7 x 10 22 s followed by a second rapid step 
up to 1.7 x 10 23 cm -2 approximately 1500 s later. After an 
additional 1500 s the absorption then rapidly decreases back 
to the base level observed in the earlier parts of the observa- 
tion whilst, almost simultaneous, the flux begins to rise into 
the flaring activity observed at the end of the observation. 

To understand this evolution we first consider the vari- 
ation in n_sr observed as the NS moves through its orbit. 
|Sidoli et al.| | |20"l2| | proposed that the NS in IGR 16418-4532 
could be accreting through Transitional Roche Lobe Over- 
flow (TRLO) due to the systems narrow orbit and the flux 
variations observed in XMM-Newton data taken in 2011. 
Such a transitional scenario was investigated by |Blondin^1 
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Stevens, & Kallman (19911 where model orbital tih pro- 



files were produced for wind accretion powered binaries with 
narrow, circular orbits (Figure 2 of that work). The overall 
shape of the profile is dominated by the line of sight through 
the bow shock about the NS with an additional, prominent 
component being generated by the tidal gas stream extend- 
ing from the LI point. The PEO was separated from the 
middle of the eclipse by an orbital phase interval of ~ 0.2 
which should, according to [Blondin, Steve ns^ fc Kallman| 
( 1991 1, result in a smooth decline in njf . This suggests that 



the increase in iih observed is not generated by a coherent 
structure in the NS environment (i.e. a tidal stream) but is 
instead more likely to be generated by an individual dense 
clump of material travelling through the line of sight to the 
NS at this time. Significant variation in the photon index 
of the powerlaw was not observed during the period of en- 
hanced uh and the unabsorbed flux remained at a constant 
level during this time (Fig. [9] red line). This is likely a sig- 
nature of absorption of a constant intrinsic flux by optically 
thin material, as opposed to an optically thick/ionised ab- 
sorber which might be expected in the case of obscuration 
by a coherent structure such as a tidal stream. We note that 
sensitive, soft X-ray observations at a later orbital phase 
would be required to further test for the presence of a tidal 
stream arising from TRLO in IGR J16418-4532. 

A dynamical estimate of the extent of the over dense 
region required to generate the excess n^ observed can 
be made by utilising the calculated relative velocity be- 
tween the NS orbital motion and the radial stellar wind (see 
above) . The increase in tih is observed for ~ 2500 s which re- 
lates to the passage of a wind clump of radius ~8x 10 10 cm 
travelling at a velocity of 470 km s -1 . This radius is consis- 
tent with the range of clump radii derived from the theo- 



retical considerations of Ducci et al. ( 2009 1 and also similar 



to the observational results of Bozzo et al. (20111 who ex- 



plained an outburst of IGR J18410— 0535 as resulting from 
the accretion of a clump of approximate radius 8 x 10 1 cm. 
Therefore we conclude that this absorption feature likely re- 
sults from an optically thin clump of stellar wind material 
obscuring the intrinsic emission from the neutron star. 



4-3.2 Absorption and flaring link? 

The temporal relationship between the tih increase and the 
onset of the flaring activity at the end of the PEO is of in- 
terest. Figure [9] shows the tih evolution together with the 
unabsorbed 0.5 — lOkeV flux derived from the fixed T fit 
to each spectrum. It can be seen that the unabsorbed flux 
only begins to rise above the base level in the observation 
as the riH returns to its pre-enhancement level. The tim- 
ing of these two events is intriguing, but on the basis of 
this dataset alone no conclusive physical connection can be 
drawn between these features and instead they may simply 
be generated by two separate, unrelated events that ran- 
domly occurred sequentially. Given the tighter constraints 
on IGR J16418— 4514s orbital geometry through the X-ray 
eclipse however, a consideration as to whether such a delay 
could be generated through physical means and indicate a 
causal link between the tih and flux evolutions observed can 
be undertaken. 

The lack of an overlap in the evolution of these param- 
eters appears to be in contradiction to the best example of 



clump accretion in an SFXT as presented in |Bozzo et al.| 
( |2011[ ). In this work an outburst of IGR J18410-0535 was 
described through the accretion of a dense clump of wind 
material. The tih was observed to rise throughout the entire 
flux evolution, only returning to its pre-flare level when the 
flux also returned to pre-flare level, suggesting that there 
was excess material in the immediate vicinity of, and line 
of sight to, the NS throughout the duration of the flare. 
However in IGR J16418— 4532 it appears that if the denser 
clump material is responsible for generating the subsequent 
flare then the clump was only in our line of sight to the neu- 
tron star prior to the flux increase. This discrepancy seems 
somewhat unintuitive but can start to be understood from 
the more detailed knowledge of the orbital geometry of the 
IGR J16418— 4532 system as a result of the X-ray eclipse. 
Using Fig. [l] we can centre the eclipse at an orbital phase 
of approximately <j> = 0.55 meaning that in the orbital phase 
range </> = 0.3 to 0.8 the neutron star is more distant than 
the supergiant companion. Hence in this phase range, which 
the PEO inhabits, stellar wind clumps that are accreted will 
only be in our line of sight to the neutron star as they are 
approaching and interacting with it. As the clump moves 
beyond the X-ray emitting region the absorption detected 
along our line of sight would then drop rapidly, as is ob- 
served in Fig. [9] as the clump moves away behind the neu- 
tron star and we are therefore unable to continue to detect 
its presence. 

To further investigate whether this geometrical effect 
could be behind the observed evolution we consider the 
free fall time of stellar wind material from the accretion 
radius. We again assume the accretion radius is given 
by R a = 2GMjvs/ u re!i an d the free fall time is given by 
tff = (R„/2GMjvs) 3 ■ Under these assumptions and using 
the previously derived system parameters, the accretion ra- 
dius is calculated as ~ 9x 10 10 cm and the associated free fall 
time is ~1500 s. This value should be considered a lower limit 
on the true infall time however as the descending material 
does not fall through the gravitational potential unimpeded, 
likely becoming a turbulent flow at some point which will 
increase the infall time. Given that the calculated free fall 
time is comparable to the observed delay between the onsets 
(and peaks) of the iih and unabsorbed flux evolutions, Fig. 
[9] and the X-ray generation region is close to the neutron 
star, it is possible that in this case the dense wind clump 
passed through the accretion sphere before the additional 
accreted material had reached the X-ray generation region, 
hence creating the delay observed in Fig. [9] 

Through considerations of the orbital geometry of IGR 
J16418— 4532, the temporal coincidence (of the order of the 
free fall time from R a ) and approximately similar durations 
of the iih and flux evolutions, we suggest that a clump of 
over dense, optically thin wind material approaching the 
neutron star at a distance comparable to the accretion ra- 
dius would be consistent with the production of both fea- 
tures. In this case IGR J16418— 4532 is showing signatures 
of accretion through the 'clumpy wind' model (in't Zand 
120051 iWalter & Zurita Heras||2007|) viewed at an approxi- 



mately edge-on inclination. However we again stress that at 
the current time it is not possible to definitively conclude 
that these two features are causally related and that they 
might instead be caused by two physically unrelated events. 
Long term monitoring of IGR J16418— 4532 across all orbital 
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phases would allow firmer conclusions to be drawn as to any 
physical link between these features through the detection, 
or otherwise, of additional uh - flux delays and their distri- 
bution in orbital phase. If it were the case that these fea- 
tures are physically related, it cannot currently be discerned 
if the difference in the spectral evolution observed in IGR 
J16418-4532 comp ared to the example of IGR J18410-0535 
i Bozzo et al.|201l"j ) is due to the same geometrical effects, as 
the orbital configuration of IGR J18410— 0535 is currently 
unknown. 

Finally there is some evidence of a slow decrease in 
the base uh level at times greater than t = 10,000 s which 
could be the result of increased X-ray photoionisation of the 
absorbing material by the enhanced X-ray flux produced 



during the flare (Hatchett & McCray 19771, however the 



observation ended before this effect could be more firmly 
identified. 



4-3.3 The nature of the 'dip' 

We now consider the properties of, and possible physical 
mechanisms behind, the X-ray dip observed towards the 
start of the PEO. Due to the unfortunate lack of EPIC-pn 
coverage of the majority of the dip region we were unable to 
perform a detailed spectral analysis of the ingress and mid- 
dip regions, as shown by the large error bars and poor x 2 
values in Figs. [7] and [8] but did achieve more sensitive cov- 
erage during the dip egress. The main observable features of 
the X-ray dip are: 

— a total duration of ~ 1500 s between t=800 and 2300 s 
that is comparable to the spin period of the neutron star 

— a sharp drop in flux lasting less than 300 s at the start 
of the feature 

— a steady, low level flux lasting for ~ 900 s with an un- 
absorbed flux of 2. Ox 10" 11 erg cm" 2 s" 1 (0.5 - lOkeV) cor- 
responding to a luminosity of 4.1 x 10 35 ergs -1 at a distance 
of 13kpc 

— a slightly less sharp return to the pre-dip flux level over 
~300 s with some evidence of the xih decreasing back to the 
pre-dip level as the flux recovers (see Fig. [9| 

— a decrease in flux observed across the full 0.2— lOkeV 
energy range during which spectral changes could not be 
identified 

— a smooth spectrum with no evidence of Fe emission 
lines. 



The presence of X-ray intensity dips, also referred to 
as 'off-states', in the light curves of wind- fed HMXB pul- 
sars is a rare phenomenon, having only been observed in 
the classical Sg- XRBs Vela X-l (|Kreykenbohm et al.|2008|), 



Centaurus X-3 (|Naik, Paul, fe Ali|2011| ), 4U 1907+ 09 (|in 't 
Zand, Baykal, fe Strohmayer|1998j) and GX 301 -2 (|Gogu§, 



Kreykenbohm, & Belloni 20111. Bozzo et al. (20121 also 



reported a similar feature in the non-pulsating, candidate 
SFXT IGR J 16328-4726 but a detailed analysis and char- 
acterisation was prevented due to statistical constraints. Dif- 
ferent mechanisms have been proposed to explain the X-ray 
dips observed in the above sources. Here we consider their 
applicability to IGR J1618— 4532, whilst also noting that 
additional, as yet unexplained, mechanisms may be acting 
to create this feature. 



The eclipsing, Roche-Lobe filling Sg-XRB system Cen 
X-3 displayed several X-ray dips during the course of a long 
Suzaku observation that covered almost an entire binary or- 
bit ( |Naik, Paul, fe Ali|2011[ ). The dips were observed up to 
energies of 40 keV and the spectra characterised by a typical 
accreting HMXB pulsar continuum with multiple Fe emis- 
sion lines that were detected throughout the observation. 
The equivalent widths (EW) of the Fe lines showed signif- 
icant enhancement during both the dip and binary eclipse 
regions of the exposure. By following the evolution of the 
Fe line EWs and the amount of absorbing material in the 
line of sight (tih values increased by a factor of 100 dur- 
ing these times), the authors concluded that the dips were 
most likely caused by the obscuration of the neutron star by 
dense structures in the outer edge of the accretion disc in 
the system. In this way Cen X-3 is a HMXB analogy of the 
'dipper' class of LMXBs. Obscuration due to dense struc- 
tures in the vicinity of the neutron star is unlikely in the 
case of IGR J16418— 4532 however. Observationally the de- 
tection of prominent Fe emission during the MEO but not 
during the dip suggests that obscuration by dense, compton 
thick material is not the cause of the intensity dip observed 
in these observations. Additionally it is not clear how such 
dense material could inhabit the IGR J16418— 4532 system 
as the supergiant is not believed to overfill its Roche-Lobe 
and therefore an accretion disc would not be present. In 
the TRLO regime such dense material could be found as 
a result of either the gravitationally focused 'tidal stream' 
coming from the LI point or as a result of a transient accre- 
tion disc that may form. However, as discussed previously, 
the orbital phase location of the dip is inconsistent with 
what would be expected if it were to originate from a tidal 



stream (Blondin, Stevens, & Kallman 19911. Sidoli et al. 
( |2012 1 suggested the presence of a transient accretion disc 
to explain the quasi-periodic flares observed in a previous 
observation of IGR J16418— 4532, however as no such fea- 
tures are present in these observations it is unlikely that a 
transient accretion disc was present to obscure the NS at this 
time. We note however that these observations cover an in- 
sufficient fraction of orbital phase to consider fully whether 
TRLO is occurring in IGR J16418-4532. 

The signature of a typical stellar wind clump obscuring 
the neutron star in IGR J16418— 4532 is believed to have 
been observed later in the PEO, as discussed previuosly, and 
displayed a very different evolution to that observed during 
the dip feature. Hence it is also unlikely that the passage of 
a single stellar wind clump through the line of sight is re- 
sponsible for the observed dip. Finally the passage of ionised 
material through the line of sight to the neutron star could 
cause a strong absorbing effect across a wide energy range 
as observed. Again however there was insufficient signal to 
test for this effect spectrally and the lack of Fe emission lines 
in the X-ray spectrum of the dip suggests this is an unlikely 
cause of the observed feature. 

The alternative explanation given for the presence of X- 
ray intensity dips or 'off-states' in the light curves of wind- 
fed SgXRBs is the action of a barrier to the accretion flow. 
|Kreykenbohm et al.| ( |2008| ) invoked centrifugal barriers to 
explain the off-states observed in Vela X-l, where the source 
fell below the detection threshold in INTEGRAL /IBIS data, 
via the onset of the supersonic propellor mechanism ( |Illar-| 
ionov & Sunyaev 1975 1 to inhibit accretion during these 
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times. More recently however Doroshenko et al. (2011 1 used 
sensitive, broadband Suzaku observations to identify the 
continued presence of X-ray pulsations during the dips of 
Vela X— 1 and argued that the dips are the result of a 
magnetic barrier that reduces the emitted X-ray flux as 
direct accretion is impeded by the neutron star magneto- 
sphere (subsonic propellor regime, see 



Bozzo et al. 20081 



ct al. (2011) used the detailed accretion scenario model of 



In this state accretion continues at a lower level through 
the Kelvin-Helmholtz Instability (KHI) which generates the 
lower flux levels detected during the X-ray dips. Doroshenko 



Bozzo et al. 



( 2008 1 to infer that a neutron star with a mag- 



netic field strength in the range of a few 10 la to ~1 x 10 14 G 
is required to produce the X-ray luminosity observed under 
this model and argue that X-ray intensity dips in wind-fed 
HMXBs are evidence of a highly magnetised neutron star 
being present in the system. The same authors also found 



si milar r esults for the HMXB 4U 1907+09 (Doroshenko et 
aTp012). 



Shakura, Postnov, & Hjalmarsdotter (2013) however, 
proposed that the X-ray dips observed in Vela X-l, GX 
301—2 and 4U 1907+09 were instead caused by a switch 
in the polar beam configuration due to a variation in 
the optical depth of the accretion column. In this model 
both states occur within a 'quasi-spherical' subsonic accre- 
tion regime (Shakura et al. 20121. The Compton cooling 



regime occurs when the source luminosity is in excess of 
~5xl0 35 ergs -1 and the optical depth of the accretion col- 
umn is high, creating a fan beam due to high levels of lat- 
eral scattering. The fan beam effectively cools plasma in 
the equatorial regions of the magnetosphere through Comp- 
ton processes, allowing material to enter the magnetosphere 
through Rayleigh- Taylor instabilities. It is this state that 
generates the 'normal' luminosity level of the lower luminos- 
ity, persistent wind-fed HMXBs (i.e. Vela X-l). The radia- 
tive cooling regime occurs when the luminosity is less than 
~5xl0 35 ergs -1 , causing a drop in the optical depth of the 
accretion column and the onset of pencil beam emission. 
In this configuration the equatorial region of the magneto- 
sphere is not illuminated and plasma can only thread onto 
field lines at a lower rate through radiative plasma cooling, 
which is independent of the level of illumination from the 
central source. The authors argue it is the transition into 
this radiatively cooled accretion state from the Compton 
cooled state that generates the observed X-ray dips. Sources 
would only move into the direct accretion scenario (super- 
sonic Bondi accretion) when the source luminosity increased 
to greater than several 10 36 ergs -1 and in-falling material is 
rapidly cooled by Compton processes before reaching the 
magnetosphere (for a full description of the quasi-spherical 



accretion model see Shakura et al. ( 2012 1 and Shakura, Post- 



nov, & Hjalmarsdotter (20131). 

Given the similarities in the observable features and lu- 
minosities of the X-ray dips observed in Vela X-l and 4U 
1907+09 to the dip observed here, we consider the applica- 



bility of the Doroshenko et al. (2011 1 and Shakura, Postnov, 



& Hjalmarsdotter ( 2013 1 interpretations of the origins of X- 
ray dips to IGR J16418-4532. We first follow the consid- 
erations of Doroshenko et al. (2011), using the theoretical 



framework of Bozzo et al. (2008), to calculate the magnetic 



field strength required to produce the observed dip luminos- 
ity (4.1x 10 35 ergs -1 ) for the range of stellar wind parame- 



ters considered in Section [4. 3| We use Equation 21 of |Bozzo] 
et al. (2008) and the conventions therein, namely: 



L KH =1.8 x lQ A5 7 lKH P a3 1 Bii IW M. 6 a^ d v s 1 



[1 + ldRaio/^RMio)] 



3/2 VPi/Pe 



(2) 



(1 + pi/pe) 



ergs 



Here the pi / p e term is the ratio of the plasma densities inside 
and outside of the magnetosphere and Rmw is the magnetic 
radius in units of 10 10 cm and defined as: 



Rm = 2 x lO'V^M 



2/7 8/7 4/7 



-6 



P33 Cm 



(3) 



where ^33 is the magnetic moment in units of 10 33 Gcm 3 . 
We calculate the magnetic field required to produce the dip 
luminosity in the cases of pi/p e = 1-0, 0.5 and 0.1 where 
a Pi/pe value of 1 gives an upper limit on the mass flow 
generated through the Kelvin-Helmholtz instability. 

Using the calculated magnetic field strengths we then 
estimate the critical accretion rate defining the transition 
between the subsonic propellor and direct accretion regimes. 
According to |Bozzo et al.| ( |2008[ ) this critical accretion rate 
is given by: 



M ilm _ 6 =2.8 x W 2 P s3 3 a 2 10d v s 

,5/2 



R uw\ 1 + l6Raio/(5R M io) 



-3/2 -1 

' gs 



(4) 



These critical values are compared to the accretion 
rates of 2.2xl0 15 gs -1 and l.lxl0 16 gs -1 (using L x = 
GMtfsMc av t/ Rns) required to generate the observed dip 
luminosity of 4.1xl0 35 ergs -1 and the pre/post-dip lumi- 
nosity of ~2x 10 36 ergs -1 respectively to evaluate if a tran- 
sition between the accretion regimes of this model could have 
occurred. Table [3] outlines the results of these calculations 
for the range of possible stellar wind velocities considered in 
Section g21 

It is seen from Table [3] that transitions between these 
two accretion regimes cannot be responsible for the observed 
dip in the light curve of IGR J16418— 4532 as the critical 
transition accretion rates are in excess of the dip, pre/post- 
dip and peak-flare (2.8 xlO 16 gs -1 ) accretion rates required 
to generate the observed luminosity levels across the whole 
parameter space investigated. Hence this suggests that a 
transition between the 'subsonic propellor' and direct accre- 
tion regimes is not causing the X-ray dip observed in IGR 
J16418-4532. 

Given the earlier arguments against obscuration as the 
source of the dip and the incompatibility of the required 
accretion rates for a subsonic propellor to direct accretion 
transition, we investigate a state change within the 'quasi- 



spherical' accretion model (Shakura, Postnov, & Hjalmars- 
|dotter||20T3"l ) as a possible cause of the X-ray dip in IGR 
J16418-4532. As IGR J16418-4532 appears to stay in a 
subsonic quasi-spherical accretion regime throughout the 
observation this suggests the X-ray dip may be generated 
through a transition between the higher luminosity fan beam 
and lower luminosity pencil beam dominated states of the 
model. Under this scenario we can estimate the magnetic 
field strength of the NS through the luminosity detected 
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Table 3. Calculated magnetic fields required to generate the observed dip luminosity in the 'Subsonic Propellor' accretion mode according 
to |Bozzo et al.j (2008). The wind condition parameter space is drawn from the limits on the terminal wind velocity range of |Romano| 



et al. 1 2012 I and the total stellar wind mass loss rate calculated using the M w /voa=9.3x 10 — Mq km parameter for each terminal 
wind velocity (see Sect, |4~3| l. The 'Transition?' columns indicate whether a transition from the subsonic propellor regime to the direct 
accretion regime is likely to have occurred for the X-ray dip and flare observed in the observation. The accretion rates required to produce 
the dip, pre/post-dip and peak flare luminosities are 2.2 xlO 15 gs -1 , 1.1 xlO 16 gs _1 and 2.8 XlO 16 gs" 1 respectively. 



pi/pe B(G) Mn m (gs 1 ) Dip Transition? Flare Transition? 

•Uoo = 800kms~\ M w =2.3xl0 -7 M yr" 1 , v ret = 571 krns" 1 

1.0 2.8xl0 13 1.2xl0 17 N N 

0.5 3.0xl0 13 1.4xl0 17 N N 

0.1 5.3xl0 13 5.2xl0 17 N N 

Voc = 1000 km s" 1 , M w =2.9xl0 -7 Mq yr" 1 , v rel = 637 km s" 1 
1.0 3.2xl0 13 3.2xl0 17 N N 

0.5 3.4xl0 13 3.8xl0 17 N N 

0.1 6.1xl0 13 1.4xl0 18 N N 

-1 



Voo = 1300 km s _1 , M w =3.8xl0 -7 M yr" 1 , v rel =746 km s~ 

1.0 4.1xl0 13 1.5xl0 18 N N 

0.5 4.3xl0 13 1.8xl0 18 N N 

0.1 7.6xl0 13 6.4xl0 18 N N 



during the dip using Equation 22 of |Shakura, Postnov, fc| 
Hjalmarsdotter (20131, namely: 



(5) 



r 1n 35 7/33 -1 

L X ,rad~W M30 er S S 

where nso is the magnetic moment in units of 10 30 
Gem 3 . The dip luminosity of 4.1 x 10 35 erg s _1 corresponds 
to a NS B-field of ~2xl0 14 G, indicating that the presence 
of a highly magnetised NS in a SFXT would be required 
by this model. Additionally the shell of material that would 
be surrounding the NS under this model could be providing 
a significant contribution to the high base level of absorp- 
tion observed in the PEO (~7x 10 22 cm" 2 ). As the derived 
B-field value is determined by the luminosity of the source 
during the X-ray dip we consider the effect of the unknown 
source distance on this interpretation. The confirmation of 
the companion as a supergiant requires the source to be at a 



large distance for all possible stellar radii ( Chaty et al. 2008 1, 



however a distance of 13 kpc is still somewhat arbitrary and 
is used here to maintain consistency with previous works in 
the literature. If IGR J16418— 4532 is located at a distance 
of 10 kpc (for example) then the dip luminosity would be 
2.4xl0 35 ergs" 1 and the implied NS B-field would have a 
value of 1.6 xlO 13 G. Whilst in this case the B-field doesn't 
reach magnetar field strengths it is nevertheless in the highly 
magnetised regime, being in excess of the nominal ~10 12 G 
value observed in many HMXB pulsars. 

Due to the lack of EPIC-pn coverage during the early 
regions of the dip it is difficult to evaluate what may have 
caused the on-set of an accretion regime transition for the 
duration of the dip. However the hydrodynamical simula- 
tions of |Oskinova, Feldmeier, fc Kretschmar| ( |2012[ ) showed 
that, whilst on average the velocity of the stellar wind fol- 
lows the velocity law stated above (Eq. [I]), strong veloc- 
ity jumps over small radial extents are also present. These 
jumps, which can have negative gradients, result from un- 
stable growth in the line-driven wind and are at their most 
prominent at separations within a few stellar radii of the 
supergiant surface. As a result of its short orbit the NS in 



IGR J16418— 4532 is occupying the most turbulent region 
of the supergiant wind and an accretion regime transition 
may have been triggered by a strong, localised variation in 
the stellar wind environment interacting with the magneto- 
sphere. 

The previous paragraphs have outlined the applicability 
of the mechanisms currently proposed to generate the X-ray 
dips observed from some SgXRBs. Through these consider- 
ations we conclude that the most applicable of these mech- 
anisms is that of an accretion regime transition between 
the two subsonic states of the 'Quasi-spherical accretion' 



model (Shakura et al. 2012b. An important implication of 



this model is that of the requirement of a highly magnetised 
neutron star to be present within the IGR J16418— 4532 sys- 
tem. Caution must be taken with these interpretations how- 
ever due to the large number of assumptions required to per- 
form these calculations. The detection and characterisation 
of additional dips is required, along with independent con- 
straints on the orbital, stellar and stellar wind parameters of 
IGR J16418— 4532, to allow the testing of the applicability of 
all of the possible causes of the dip to the IGR J16418— 4532 
system under more stringent physical constraints. At the 
current time it is therefore not possible to definitively iden- 
tify the physical mechanisms that create the dip feature. 



4.4 The nature of IGR J16418 4532 

The strongly variable stellar wind environment being expe- 
rienced by the NS in IGR J16418— 4532 may also explain 
the origin of the enhanced dynamic range (> 10 2 ) observed 
in this system compared to classical Sg-XRBs. Many Sg- 
XRBs have longer orbits and should therefore occupy less 
turbulent regions of their supergiants atmosphere resulting 
in a lower level of variability. The sharp velocity and den- 
sity variations near to the surface of the supergiant in IGR 
J16418— 4532 however may be helping to generate the ob- 
served SFXT level of variability in the system through their 
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influence on the size of the accretion radius long with the 
achieved mass capture rate. 

The transient nature of IGR J16418-4532 is defined 
more by the large source distance of ~13 kpc however. Figure 
[I] shows that outside of the eclipse IGR J16418— 4532 is de- 
tected as a weak persistent source in the INTEGRAL /IBIS 
18— 60keV band (where most SFXTs were discovered) at 
a count rate of 1.2 counts per second (~7mCrab). If IGR 
J16418— 4532 was located at a distance of 2 kpc however 



(comparable with Vela X-l, Nagase 19891 it would be de- 



tected as a bright persistent Sg-XRB source with an unusu- 
ally high dynamic range and an average flux of ~300 mCrab 
in the INTEGRAL band. Hence it would not be viewed as a 
transient source. Therefore one may take the view that IGR 
J16418— 4532 is a classical Sg-XRB with an enhanced X-ray 
dynamic range resulting from the turbulent stellar wind en- 
vironment being experienced by the NS as a result of the 
short orbit orbital period of the system. Its transient na- 
ture is then being generated as a result of its large distance 
making the persistent emission undetectable in single IN- 
TEGRAL/IBIS or Swift /BAT observations, but where the 
flares are instead detected as transient outbursts. This is at 
odds with many other SFXTs where the driver of the ex- 
treme level of flux variation is the large scale variation in 
stellar wind environment experienced by a NS over a longer, 



more eccentric orbit (e.g. SAX J1818.6-1703, Bird et al 



20091. Shakura et al. (20121 also suggest that the difference 



between SFXTs and classical wind-fed Sg-XRBs could be 
understood through the proportion of time spent in the fan 
beam dominated (higher luminosity) and pencil beam domi- 
nated (lower luminosity) regimes within the subsonic, quasi- 
spherical accretion mode. Classical Sg-XRBs (e.g. Vela X-l) 
spend the majority of the time in the fan beam dominated 
regime and occasionally enter the pencil beam regime dur- 
ing X-ray dips. Whereas SFXTs spend the majority of their 
time in the pencil beam dominated regimes and occasion- 
ally enter the fan beam regime during short periods of in- 
creased accretion, i.e. outbursts. Given IGR J16418— 4532s 
emission history and the observations presented here this 
further suggests that, fundamentally, IGR J16418— 4532 is a 
classical wind-fed Sg-XRB displaying an enhanced dynamic 
range and is observed as a transient source only due to its 
large distance. 

This difference in the fundamental cause of the variabil- 
ity observed in different SFXT systems is suggestive of the 
fact that the SFXT class is a phenomenological group de- 
fined by the duration and dynamic range of the flares and 
consists of multiple populations of wind-fed Sg-XRBs whose 
enhanced level of variability arises from different physical 
origins. These populations could also be the result of dif- 
ferent evolutionary paths followed to reach the current Sg- 



XRB state as Liu, Chaty, & Yan (20111 argued that some 



longer orbital period SFXTs (namely IGR J18418-0311 and 
IGR J11215— 5952) were most likely to have evolved from 
OeXRB systems. Such systems, as with the BeXRBs known 
currently, would likely have had longer, more eccentric or- 
bits which would not change significantly when the main 
sequence Oe star evolved off the main sequence, thus provid- 
ing the current orbital configuration inferred in the SFXTs 
( Liu, Chaty, fc Yan|201l] ). 

To gain further understanding into the accretion pro- 
cesses occurring and the possible different sub-populations 



that make up the SFXT class will require high sensitivity 
X-ray observations covering wide regions of orbital phase. 
Additionally an accurate characterisation of the supergiant 
stellar and stellar wind parameters in individual SFXTs will 
aid in identifying the accretion scenarios occurring under 
stricter physical constraints. In particular achieving sensitive 
monitoring of SFXTs over the entire orbital period range 
is imperative in identifying if there is a systematic change 
in behaviour between sources whose orbital parameters are 
more consistent with Sg-XRBs and those whose parameters 
more closely resemble BeXRBs. 



5 CONCLUSIONS 

In this work we have presented new combined INTEGRAL 
and XMM-Newton observations of the intermediate SFXT 
IGR J16418-4532. We have identified features in the XMM- 
Newton data that are consistent with the presence of dense 
stellar wind clumps as predicted by the nominal 'clumpy- 
wind' model of SFXT emission along with the first ever de- 
tection of an X-ray 'dip' in the light curve of a pulsating 
SFXT. Through considerations of the cause of this dip and 
the luminosity levels detected we conclude that the most 
likely cause was the switching of accretion modes from a 
Compton cooling dominated to a radiative cooling dom- 
inated regime within the framework of a subsonic quasi- 
spherical accretion flow (Sha kura et al.||2012[ ). Further ob- 
servations and characterisation of IGR J16418— 4532 are re- 
quired to definitively identify this scenario as the correct ex- 
planation of the observed behaviour but, if confirmed, IGR 
J16418— 4532 would be required to host a highly magne- 
tised, B ~ 10 14 G, neutron star. In this case the observa- 
tions presented in this paper would represent the first ob- 
servational evidence of such an object being hosted in an 
SFXT. It is likely that the observational properties of IGR 
J16418— 4532, that of a persistent wind-fed SgXRB with an 
enhanced dynamic range, are in part defined by the high 
level of turbulence in the stellar wind being experienced 
by the neutron star in its short orbit close to the compan- 
ion supergiant star. The orbital geometry and flare gener- 
ation mechanisms observed in IGR J16418— 4532 represent 
a departure from the nominal model of SFXT behaviour. 
This suggests that there may be sub-populations within the 
SFXT class that, whilst having some similar defining obser- 
vational characteristics (dynamic range, flare luminosities, 
stellar companions) , generate the observed X-ray behaviour 
through a variety of different processes. 
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